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The eigenvalue density of a Wilson loop matrix W associated with a 
simple loop in two-dimensional Euclidean SU(iV) Yang-Mills theory under- 
goes a phase transition at a critical size in the infinite- iV limit. The averages 
of det(z-H^)-i and det{l + uW) / {1 - vW) at finite N lead to two different 
smoothed out expressions. It is shown by a saddle-point analysis that both 
functions tend to the known singular result at infinite N. 

1. Introduction 

In two Euclidean dimensions the eigenvalue distribution of the SU(A^) 
Wilson matrix associated with a non-selfintersecting loop undergoes a phase 
transition in the infinite- limit as the loop is dilated [Ij . This phase tran- 
sition has universal properties shared across dimensions and across analog 
two-dimensional models [21 [3]. Thus, a detailed understanding of the tran- 
sition region in 2D is of relevance to crossovers from weakly to strongly 
interacting regimes in a wide class of models based on doubly indexed dy- 
namical variables with symmetry SU(A^). Building upon previous work 
in El E]) new results in this context have been obtained in [7]. Some of 
these results are presented here. 

We are focusing on the eigenvalues of the Wilson loop. The asso- 
ciated observables are two different functions of an an- 
gular variable 6. At infinite N the two functions have identical limits: 

For a specific critical scale, the nonnegative function PooiO) exhibits a 
transition at which a gap centered at 9 = ztyr, present for small loops, just 
closes. This transition was discovered by Durhuus and Olesen in 1981 [T]. 
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2. Eigenvalue densities 

The probability density for tlie Wilson loop matrix W is given by the 
heat kernel (see for example [8j and original references therein) 

VN{W,t) = J]d,Xr(W^)e-^^2(r) (1) 
r 

with t = XA, where A = g'^N is the standard 't Hooft coupling and A 
denotes the area enclosed by the loop. The sum over r is over all distinct 
irreducible representations of SU(A^) with dr denoting the dimension of r 
and C2(r) denoting the value of the quadratic Casimir on r. Xr{W) is the 
character of W in the representation r and is normalized by Xr(l) = dr- 
Averages over W at fixed t are given by 

{0{W)) = j dWVN{W,t)0{W) , (2) 

where dW is the Haar measure on SU(A^) normalized by / dW = 1. Note 
that we have / dWVN(W,t) = 1. Any class function can be averaged when 
expanded in characters using character orthogonality. 

Because in the sum over r in ([T]) each representation is accompanied by 
its complex conjugate representation, it is easy to see that 

{0{W)) = {OiW^)) = {OiW*)) , (3) 

implying identities relating (det(2; — W)), (det(z — W)~^), and (det(l + 
uW)/{l — vW)) to the same objects with z ^ 1/z, z ^ z* , u,v ^ 1/u, 1/v, 
and u,v — > u*,v*, respectively. 

The density functions p^™ and p^"^ are obtained from 



Gr(^) = >°S<d«' - ') (5) 

through {i = true, sym) 



p%{e) = 2 hm Re e^^+^G^(e'^+^) - 1 . (6) 

e— »0+ L J 

The are real on the unit circle parametrized by the angle \9\ < vr, even 
under 9 —9, and depend on the size of the loop. Both are positive 
distributions in 9, normalized by 



(7) 



3 




0.5 1.0 1.5 2.0 2.5 3.0 0.5 1.0 1.5 2.0 2.5 3.0 



Fig.l. Plots of p^"'(6',T) for T = 2 (left), T = 5 (right), and N = 
3, 5, 10, 25, 50, 100, 250 together with p^{e, T). 



Only has a natural interpretation at finite N ^ it literally is the 

eigenvalue density. If the eigenvalues of W are e*"-' with j = 0,l,...,A^ — 1, 
it is given by [7] 

pT{0) = %Y.^52A0-a,{W))) = ^(TV52.(^ + ilog(M^))>. (8) 

3 

p^^^ determines (Tr/(H^)) for any function /. 

The density /O^™ is determined by the averages of the characters of W in 
all totally symmetric representations of SU(A^). This function is of interest 
mainly because it obeys simple partial integro/differential equations which 
are exactly integrable f5] . 

/O^™ has an explicit form in terms of rapidly converging infinite sums 
[5] , which can be evaluated numerically for arbitrary N to any desired pre- 
cision. In Fig. [l] we show how p'^^iO) approaches the infinite-A'" result 
Poo{d) of Durhuus and Olesen p^™ is monotonic on each of the seg- 
ments (— 7r,0) and (0,7r) with the maximum at ^ = and the minimum at 
9 = ivr. In addition to these numerical results, it would be useful to com- 
pute analytically the asymptotic expansion of p^™(^) in (cf. section 

i- 

It turns out that the appropriate area variable for is not t = XA 
but r = t(l - l/AT) ig . When T) is compared to p%''%0, t) , the 1 /N 

correction in t relative to T has to be taken into account (so far, the size 
dependence of the has been suppressed) . 

The infinite-A^ critical point is at T = 4. For T > 4, p^™(6l,r) ap- 
proaches poo{0,T) by power corrections in |5, . For T < 4, poo{0,T) is 
zero for \e\ > OdT), where < OdT) < vr and 6*^(4) = vr (cf. Fig.[TJ. In this 
interval /0^™(6', T) approaches zero by corrections that are exponentially 
suppressed in N. 
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Fig. 2. Plots of the densities p%'"'{e,t) (red, solid) and p%"'{9,T) (blue, dashed) 
for t = 2 (left) and t = 5 (right), = 10 (top), and = 50 (bottom). 



The true eigenvalue density p^™'^{6,t) has N peaks (in the interval 
[-7r,7r]) and oscillates around p'^ {e,T) (cf. ). Plots of p%'"'{e,t) and 
p^™(0, T) are shown in Fig. [2] for t = 2 and t = 5. 

3. Integral representations 

The density p^™^ can be obtained from the expectation value of [7] 

p=0 (}=0 

with |f I < 1. Xp{W) (resp. Xqi^)) denotes the character of in a totally 
antisymmetric (resp. symmetric) representation. When we set u = —v + e 
and expand to linear order in e, the LHS reads 

R{-v + e,v,W) = l-eTT-^^. (10) 

After decomposing the tensor product ® into irreducible representa- 
tions, we obtain for the expectation value of the trace (due to character 
orthogonality) 

^ ' p=0 0=0 
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where C{p, q) and d{p, q) denote the value of the quadratic Casimir and the 
dimension of the irreducible representation identified by the Young diagram 





1 


2 






q 


1 
p 





(12) 



and are given by [9] 



C{p, q) = {p + q+l){N- ^^l'^^ +q-p 



1 



N 



d^p) 



d'iq) 



N p+q+l 
N + q-V 

q 



(13) 
(14) 
(15) 



We can exactly calculate sums of the form (with \z\ < 1) 



N-l 



Y,uPd^{p){N -p) = N{l+u)^-\ ^z''d^{q){N + . 

p=0 q=0 



N 



(1_ 2)^+1 • 
(16) 



To factorize the sums over p and g in (11 ), we first write 

1 



dpff 



(17) 



p + q + l 

The t-dependent weight factor is the exponent of a bilinear form in p and q 
(given by ( 13 )) . By a Hubbard-Stratonovich transformation the dependence 
on p and q can be made linear. Performing the (independent) sums over p, 
q then leads to [7] 



R{v) 



e 2 

t 

[l — vpe 



dxdy 
27r 



-x-t/2 



N-l 



(18) 



(valid for \v\ < 1). Now the entire dependence on A'' is explicit. The infinite- 
N limit of /O^"*^ can be obtained from this integral representation by using 
a saddle-point approximation for the integrals over x and y (cf. Sec. [5]). 
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4. Asymptotic expansion of 

{det{z - Wy 



An integral representation for ^'•^•'(z) 
aes /O^™, is obtained in e 
integral is needed (valid for |2:| > 1 



which deter- 



mines /O^™, is obtained in a similar manner [5]. In this case, only a single 



N 
2^ 



du e 



2T " 



ze 



-N 



(19) 



The aim of this section is to construct an asymptotic expansion of /3^™(0, T) 
in powers of 1/A (cf. [3). To this end we perform a saddle-point analysis 
of the integral in (19), from which P/^™ can be obtained via ^ and ([6]). 



.1. Saddle-point analysis 



1 the integrand of (19) has singularities on the real-n axis. We 
,e+ie ^j^gj-g e > ensures that Izl > 1 but will later be 



For \z 
therefore set z 

taken to zero. The integrand of (19) can be written as exp(— A/(n)) with 



u 

2f 



+ log 



ze 



6 2 



(20) 



We now look for saddle points of the integrand in the complex-u plane, 
which we label by u = iTU{e, T), where U{e, T) = Ur{9, T) + iUi{e, T) is a 
complex-valued function of 6 and T. The saddle-point equation turns out 
to be 



t+ie 



(21) 



-TU(e.T) U{0^T) + l/2 
U{9,T) - 1/2 

For e = 0, this is equation (5.49) in [lOj and is related to the inviscid complex 
Burgers equation via equation (5.44) there. In the present notation, the 
latter equation has the form 



dU .^BU 



0. 



Taking the absolute value of (21 ) leads to the equation 
C/2 = Ur co\h{TUr + e)-U'^ -\. 



(22) 



(23) 



For e = 0, this equation has been investigated previously in ^O]- However, 
here we keep e > for the time being. The singularities of the integrand 



of (19) then all have Ur < 0. Equation (23) describes one or more curves 



in the complex- [/ plane on which the saddle points have to lie (for a given 
value of 9, the saddles are isolated points on these curves). 



7 




o 



o 



Fig. 3. Examples of the contours in the complex-?/ plane described by equation 



([23| for r = 3 (left), T = 4 (middle), and T = 5 (right). The red dashed curves 
are for small e > 0, while the black solid curves are for e = 0. For our saddle-point 
analysis we keep e > 0. 



In Fig. [3] we show typical examples for these curves for T < 4, T = 4, 
and r > 4, where e has been chosen sufficiently close to zero. (The closed 
contours always enclose the points U = 1/2 or U = —1/2. For T > 4 and 
larger e, the closed contour in the left half-plane would be missing, but right 
now we are not concerned with this since we are only interested in the limit 



e 0^.) Analyzing (21) numerically we find, for all values of T, that for a 
given value of 9 there is always one (and only one) saddle point on the closed 
contour in the right half- plane, i.e., with Ur > 0. Note that we are showing 
the complex-C/ plane, in which the original integration contour corresponds 
to the imaginary axis. The integration contour can be smoothly deformed 
to go through the (single) saddle point in the right half-plane along a path 
of steepest descent. No singularities are crossed since they all have Ur < 0. 
There are also saddle points on the contour(s) in the left half-plane (in 
fact, there are infinitely many on the open contour), but these need not be 
considered. Figure |4] shows an example for the location of the saddle points 
and the deformation of the integration contour in the complex- ti plane. 

Once the integration contour has been deformed to go through the saddle 
point, we can safely take the limit e ^ 0"*". Parametrizing the contour 
in the vicinity of the saddle point by u = n -|- xe*^, where x is the new 
integration variable corresponding to the fluctuations around the saddle 
and fi is the angle which the path of steepest descent makes with the real-n 
axis, 'il)^^\e''^ ,T) is given, up to exponentially small corrections in A^, by 



9{x) = ^ {.e^^ + ^TUiO, T)f + logsinh - + TU [0 ,T) 

(25) 



We can now expand g{x) in x. The linear order vanishes by construction. 
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Fig. 4. Example for the location of the 
saddle points and the deformation of the 
integration contour in the complex-u plane 
for T = 5 and 6* = 3 (with small e > 0). 
The thin solid lines are lines of constant 
Re/(w), the arrows point in the direction 
of increasing Re f{u) . On each of the closed 
orange curves there is one saddle point (red 
dot and blue dot), and on the open or- 
ange curve there are infinitely many saddle 
points, but only one of them in the region 
shown in the plot (green dot). The thick 
blue curve is the integration path along the 
direction of steepest descent through the 
relevant saddle point. 



The second order gives a Gaussian integral over x, resulting in 

«W(e- T) . ,^.^I^[i;Mz£<mZ . (26) 

v^i-r(i/4-c/2(e,T)) 

Note that the factor e~*^ cannot be pulled out of the term in square brackets 
because periodicity in 6 would be lost. 

There is a potential complication. In principle, ^"(0) and therefore the 



denominator in (26) could be zero, which would mean that the integral 
over X cannot be performed in Gaussian approximation. For T > 4, it is 
straightforward to show that ^"(O) is never zero. For T < 4, one can use 



(21 ) to show that ^"(O) = only for the saddle points corresponding to the 
two angles 9 = ±0c{T) at which poo{0,T) becomes zero (see section[2]). This 
means that for \6\ = 9c(T) the asymptotic expansion in diverges, and 
that it converges ever more slowly as \9\ 9c from below. 

Note that for T < 4 and 9c(T) < l^l < vr the function p%^{9, T) is expo- 
nentially suppressed in A^. The study of the large- asymptotic behavior 
in this region requires more work. 

4-2. Leading- order result 
Equation (26) is the leading order in the expansion of ■0(^)(e*^, T). 



We now show that it leads to {9,T) poo{9,T) as A*" — > oo. We first 



write ( 26 ) in the form 



llogV(^)(e*^^) = I - /(n) + Oil/N) . (27) 
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Note that in this order we do not need the denominator in (26), which 
corresponds to f"{u) (or g"{0)). Using u = iTU this leads to 



GT{z) 



^_\_TV + ^( VA^) = i ( + ^ ) + ^( VA^) (28) 



where in the last step we have used the saddle-point equation (21 ). Equation 
(|6| then gives 



lim pl'^ie.T) =2ReU{e,T), 

N^oo 



(29) 



which equals Poci^^^T) of DO [HIIJ since U{6,T) satisfies (21) (which leads 
to (III) below with A = C/ - 1/2 and v = 1/z). 



4-3. 1/N correction to poo 

Higher-order terms in the 1/A^ expansion of i/^(^)(e'^, T) can be obtained 
in the standard way by considering higher powers of x in the expansion of 
g{x), resulting in integrals of the type J^^dx x'^"'e~^ ^^^^ with n G N. 
However, if we are only interested in the 1/A^ correction to poo{0,T) the 
result (26) is already sufficient (1/A^ corrections to this result would give 
corrections to poo{0,T)). Therefore we now write 



llog^W(e^^^) 



fin) 



2N 



log[T/"(u)] + 0(l/iV2) , (30) 



which leads to 

pT{0,T) = 2Re 



Uil + 



1 r(i/4-c/2) 
iv[i-r(i/4-c/2)]2 



+ 0{1/N'^). (31) 



Note that for T < 4 and \0\ Gc(T) (from below) the denominator of the 
1/N term approaches zero, which corresponds to the complication discussed 
[47T| Note also that for T < 4 and \e\ > 9^ the saddle point U{9, T) 



m section 



is purely imaginary so that both the leading order and the 1/A^ term are 
zero. This confirms that the above saddle-point analysis is not the right 
tool to compute finite- effects in this region. 

In Fig. [s] we show examples for the 1/A^ corrections to pao{9,T) for 
iV = 10 and r = 2 and 5. 



5. Saddle-point analysis for p]^' 



true 



We now turn to the integral representation (18) to take the first steps 
in a 1/A^ expansion of p^^^'^(9,t) (cf. [7 ). Since this integral representation 
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Fig. 5. Examples for the 1/N corrections to Poo{0,T) for TV = 10, T = 2 (left), 
and T = 5 (right). Shown are the exact result for p^^'°^{6,T) (blue dashed curve), 
the infinite- result p^{d,T) (red dotted curve), and the asymptotic expansion of 
p^^™{0,T) up to order 0{1/N) (black solid curve). We observe that the asymptotic 
expansion converges rapidly for small \6\ and more slowly for larger \9\. 



was derived for |f | < 1, we set v = e^^ 
limit e ^ at the end. We write (18) as 

Ar2 



with 1^1 < vr, e > 0, and take the 



R{v) 



t 



-e 2 



dxdy 
27r 



_ iV 

dpe 2t 



(x2+y2) + ^(x+ij/)2-i(x-iy) 







X e 



(Af-l)log(l-s;pe-^-*/2)-(Ar+l)log(l-i'pe'!'-*/2) 



(32) 



At large N , the integrals over x and y decouple at leading order and can be 
done independently by saddle-point approximations. Let us start with the 
integral over y since it is conceptually simpler. The y-dependent coefficient 
of the term in the exponent in equation (32) that is proportional to — A'^ is 



f{y) 



1 

2t 



y +log 



1 — vpe 



ly- 



(33) 



Substituting y = u — it/2 = it{U — 1/2) (with u = itU in analogy to 
section |4]) results in exactly the same integrand that was already considered 
in section |4] with the replacements T — > t and z —>■ l/vp (with \vp\ < 1) 
and with an integration over u that is now along the line from — oo + it /2 
to +00 + it/2. Since there are no singularities between this line and the 
real-u axis we can change the integration path to be along the real-n (or 
imaginary-f/) axis. Now everything goes through as in section [4| The 
saddle-point equation reads 



-tu 



U + 1/2 
U -1/2 



1 

vp 



(34) 



which is equivalent to (21). In Fig. [6] we show the contours in the complex- 
U plane on which the solutions of the saddle-point equation have to lie. 
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Fig. 6. Contours of solutions of equation (34 1 in the complex- ?7 plane at t = 3 
(left), t = 4 (middle), and t = 5 (right) for p — I (black, solid), p — 0.9 (red, 
dashed), p — 0.6 (green, dot-dashed), and p — 0.3 (blue, dotted). In the figures 
(but not in the analysis) we have taken \v\ = 1 for simplicity. 



(For sufficiently small p we now encounter the case mentioned in section 4.1 
where for t > 4 the closed contour in the left half-plane is missing.) The 
relevant saddle point, which we denote by yo{6,t,p), is again on the closed 
contour in the right half-plane. For decreasing p this contour contracts, but 
this makes no difference to our analysis. The result for the y-integral is 



given by an expression similar to (26). 



We now turn to the integral over x. The x-dependent coefficient of the 
term in the exponent in equation (32) that is proportional to — is 



fix) 



1 

2t 



log 



1 — vpe 



-t/2 



-f{ix) 



(35) 



Substituting x = —iu — t/2 = t{U — 1/2) (with u = itU) again leads to the 
integral considered in section [4] and the saddle-point equation ( 34 ) , except 
that the integration is now along the real-C/ axis. The positions of the saddle 
points of the x-integral are obtained by rotating the saddles of the y-integral 
by — 7r/2 in the complex-C/ plane, i.e., Xg = —iys- At a saddle point we have 



fix,) 



1 

t t 



1 + 



X. 



f"{ys) , 



(36) 



and therefore the directions of steepest descent through a saddle ys and the 
corresponding saddle Xs = —iys are identical (no rotation). By analyzing 
the directions along which the phase of the integrand is constant, we find 
that the integration contour can always be deformed to go through the (sin- 
gle) saddle point in the right half- plane in the direction of steepest descent. 
Depending on the parameters p, v, and t, there is either one or no additional 
saddle point on the contour(s) in the left half-plane through which we can 
also go in the direction of steepest descent. If there is such an additional 
saddle point, we find that its contribution to the integral is always exponen- 
tially suppressed in N compared to the saddle point in the right half-plane 
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Fig. 7. Example for the location of the saddle points and the deformation of the 
integration path in the complex-f7 plane for t — 5 and p — 0.95. The dashed black 
curves (two closed, one open) are the curves on which all saddle points have to lie, 
cf. (23 1. In this example 9 — 3.0. On each of the closed curves there is one saddle 



point (red dot and blue dot), and on the open curve there are infinitely many 
saddle points, but only one of them in the region shown in the plot (green dot). 
The thin solid lines are lines of constant Re f{x) and Re f{y). The arrows point in 
the direction of increasing Re f{x) or decreasing Re /(y). The dashed blue curve is 
the integration path for the y-integral along the direction of steepest descent. The 
solid red-blue curve is the integration path for the a;-integral along the direction of 
steepest descent. 



and can therefore be dropped from the saddle-point analysis. In addition, 
there are infinitely many more saddle points on the open contour in the left 
half-plane. However, we cannot deform the integration path to go through 
these points in the direction of steepest descent and therefore do not need 
to include them. An example for the location of the saddle points and the 
deformation of the integration path is given in Fig. [7| To summarize, the x- 
integral can be approximated by the contribution of the single saddle point 



in the right half- plane, which again leads to an expression similar to (26|. 

Combining the saddle-point approximations for the integrals over x and 
y, we find that, up to exponentially small corrections in A^, the integral in 



equation (32) is given by 



R{v) = e-*/2 f'dp^i -J^) ^— - e-o , (37) 

t Jo ^27r ViV/"(xo)y (l-t'pe-o-t/2)2 ^ ) 



where xq = xo{0,t,p) is the dominating saddle point of the x-integral. xq 
is a solution of the saddle-point equation obtained by differentiating f{x), 
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which can be written as 



e-.o-i/2 = (38) 



and leads to 



With (36) we obtain 



nxo)(l-^.e-o-|)^ = ^+;°(^+;o) (40) 
V / {t + xoY 

and 

R{v) = -^e-i /' dp , e-- . (41) 

Differentiating equation (psl) with respect to p leads to 



dxo ^ 1 xo{t + xq) ^ ^^-^a-tn + ^o)^ (42) 
dp /? t + (t + xo) t + 2;o (t + xo) ' 

which yields 

K^) = -T [' dp^ = ~[xo{e,t,p = l)-xo{e,t,p = 0)] . (43) 
tv Jo c*/) 

We know from (38) that XQ{6,t,p = 0) = 0. If we parametrize xo{9,t,p = 
1) = X{9, t)t, where X{9, t) has to solve 

^ = ie^(A+i/2) _ 1 ' (44) 
and take the limit e ^ 0"^, we end up with 

i?(^) = -^^MM, ^ = e^o_ (45) 

V 

Here we need to keep in mind that we have to pick the solution of equa- 
tion (44) which corresponds to the dominating saddle point xq of the x- 



integral for \vp\ < 1. 
Using (|6| and 



G(.) = ^-^i?(z-i) (46) 



we obtain 



lim pT^O, t) = 1 + 2 Re A(0, t) , (47) 

which is equal to PooiO, t) [H [TT]. Keeping higher orders in the saddle-point 
approximation (as explained in section [Xs] ), we can compute the asymptotic 
expansion of /9^"''(0,t) in powers of 
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6. Conclusions 

The probability distribution of Wilson loops in S\J{N) YM in two Eu- 
clidean dimensions can be written as a sum over irreducible representations 
(where only dimension, second-order Casimir, and character of W enter). 
This allows for the derivation of integral representations for different den- 
sity functions (including the true eigenvalue density), which have the same 
infinite- A?^ limit. These integral representations, where N enters only as a 
parameter, are exact for any finite A^. Results at infinite A'" can be obtained 
by saddle-point approximations. Next-order terms in an expansion in 
give reasonable results in the interval where Poo(^) > 0. More work is needed 
to get finite N effects where Poo{G) = and corrections are exponentially 
suppressed in N. 
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